Co-Evolution of Central Black Holes and Galaxies 

Proceedings lAU Symposium No. 267, 2010 © 2010 International Astronomical Union 

B.M. Peterson, R.S. Somerville, & T. Storchi-Bergmann, eds. DOI: 00.0000/XOOOOOOOOOOOOOOOX 

A Complete Census of AGN and Their 
Hosts from Optical Surveys? 

Vivienne Wild^, Timothy Heckman^, Paule Sonnentrucker^, 
Brent Groves^, Lee Armus'^, David Schiminovich^, 
Benjamin Johnson®, Lucimara Martins^, and Stephanie LaMassa^ 

^ Institut d'Astrophysique de Paris, 75014 Paris, France. 
Email: wild@iap.fr 

^ Dept. of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 
^ Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
* Caltech, Spitzer Science Center, MS 314-6, Pasadena, CA 91125, USA 
^ Dept. of Astronomy, Columbia University, New York, NY10027, USA 
® Institute of Astronomy, Madingley Road, Cambridge, CB3 OHA, UK 
^ NAT-Universidade Cruzeirodo Sul, Sao Paulo, Brazil 

Keywords, infrared: galaxies, galaxies: active, galaxies: starburst, galaxies: ISM 



1. Introduction 

In the last few decades, large optical spectroscopic surveys have allowed enormous progress 
to be made in studying the galaxy population, both the overall demographics and in the iden- 
tification of rare populations. The role of AGN on the evolution of the galaxy population has 
also become an increasingly common topic of research, largely driven by theoretical models and 
simulations which suggest some additional input of energy into the ISM and IGM is required in 
order to move galaxies from the blue to the red sequence and prevent star formation in massive 
galaxies. 

At low to intermediate redshifts, the sheer statistics of current optical surveys allow us to 
identify and study significant numbers of rare objects which represent transient phases in the 
lifetime of galaxies: starburst, post-starburst and "green-valley" galaxies, major mergers, and 
powerful AGN. These brief periods of major disruption in the life of a galaxy may be insignificant 
compared to a typical galaxy's total life span, but hugely significant globally, i.e., when the effects 
are integrated over the whole population and the age of the Universe. For example, summing the 
growth of black holes attained during rare, high growth-rate phases (high L/I/Edd) gives a black 
hole mass density close to that observed in the local Uni verse, implying that lower growth-rate 
event s are relatively unimportant for building black holes (|Yu fc Tremainell2002l : lHeckman et al.l 
\2004 ). 

Optical surveys allow us to locate sufficient samples of galaxies experiencing transient phases 
of evolution. There is an obvious benefit in using the information provided in the detection 
survey to study the galaxies in greater detail, but, with the complex physical processes at play 
during these periods of disruption, it is uncertain how far optical studies alone can help us in 
understanding the true impact of these transient phases on the life-cycle of the galaxy. How much 
physical information can we really infer from the optical when faced with extreme star formation, 
dust production, AGN activity and complicated spatial geometries of disturbed galaxies? 

Here, we present ongoing work comparing star formation and AGN properties derived from 
mid-infra-red spectra from the Spitzer space telescope, to those derived from optical spectra. 
We combine information from four surveys, three of which probe an "unusual" phase in the life 
of a galaxy and one control sample of "ordinary" galaxies. 
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2. Samples 

SSGSS: The cont rol sample is the Spitzer-SDSS-GALEX Spectroscopic Survey (SSGSS, see 
lO'Dowd et al.|[2009l ). a representative sample of 100 local galaxies selected from the Sloan Digital 
Sky Survey (SDSS) with GALEX photometry for follow-up spectroscopy with the Spitzer-YRS>. 
Their mass, color, star formation rate and redshift distribution are representative of galaxies in 
the SDSS survey. 

Dusty B aimer- strong: These galaxies were selected from the SDSS-DR4 to be candidate "post- 
starburst" AGN with strong Balmer absorption lines. But they also possess very high dust 
contents which may mean the str ong Balmer abso rption lines are caused by dust geometry 
rather than star formation history (jPoggianti &: Wu| [2000). 

ULIRGs : These are the 10 ULIRGs in the IRAS Bright Galaxy Sample; the Spi tzer IRS low- 
and h igh-resolution observations of these galaxies have been presented in detail in lArmus et al.l 
Of the 10 BGS ULIRGs, six are found in the SDSS spectroscopic sample, three do not 
lie within the SDSS footprint, and one is too faint to have been targeted. 

Seyferts: The narrow-line Seyfert sample is composed of the top 20 [O in] A5007 flux emitters 
out of all SDSS-DR4 main-sample galaxies classified as Seyfert galaxies using the diagnostic line 
ratio plot of [Nn]/Ha vs. [Oin]/H^. 




Figure 1. Left: The optical BPT diagram of nebular emission-line ratios. Only those galaxies 
with signal-to-noise ratio S/N >3 in all four emission lines are plotted. Right: The Spoon et 
al. (2007) IR diagnostic diagram. Only those galaxies with mean per-pixel S/N >4 in the 
wavelength range 5.4-5.9 fim axe plotted. Formal upper limits on 1^6.2 are not calculated: 14^6.2 
values < 0.05 should be considered consistent with zero. 

In Figure[T] we place the four samples into context, for both optical and IR astronomers. In the 
left-hand panel, we show the BPT diagram of optical emission-line ratios which separates galaxies 
with nebular emission dominated by star formation (left-hand branch) from those dominated by 
higher ionization processes occurring in the galactic nucleus (right-hand branch). The simplest 
explanation for the population that lies between the two demarcation lines is that both star 
formation and an AGN are present ("composite- AGN"). This is supported by their younger mean 
stellar ages as meas ured from the stellar continuum, compared to galaxies with "pure- AGN" 
emission line ratios (jKauffmann et al] l2003l : iKewlev et al.ll2006f ). The gray contours indicate 
the joint distribution of all SDSS-DR7 galaxies, the symbols show the four samples described 
above. Both the dusty-Balmer strong galaxies and the Seyferts were pre-selected to lie in the 
AGN branch. All six ULIRGs, selected purely on IR luminosity, show optical emission-line ratios 
which indicate the presence of an obscured AGN. In t he right-hand panel of Figure[T]we place all 
our samples onto the diagram of ISpoon et al] |20o3) which compares the rest frame equivalent 
width of the 6.2 ^m PAH feature caused by on-going star formation, with the strength of the 
silicate absorption feature at 10 /im caused by dust absorbing the mid-IR continuum along the 
line of sight (see following section). As expected from their optical classification as AGN, both the 
Seyfert sample and the SSGSS composite- AGN have on average lower Wq.2 than ordinary star- 
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forming galaxies, indicative of the destruction of small grain PAHs in the strong radiation field 
of the QSO. The Balmer-strong AGN show stronger silicate absorption than ordinary galaxies, 
once again qualitatively consistent with their large Balmer decrements (Ha to H/3 emission-line 
ratios) in the optical. 



3. Dust 

The effect of dust grains intercepting stellar and nuclear light complicates the analysis of 
integrated galaxy spectra. In the optical, extinction of the light along the line-of-sight can be 
quantified through a measured Balmer decrement (Ha to H/3 emission line ratio), combined with 
an empirically measured dust extinction curve (i.e. the Milky Way, LMC or SMC). However, the 
additional scattered component that makes up a full attenuation curve important for studies 
of integr ated galajcy light ha s only been measured empirically for the continua of starburst 
galaxies (Calz etti et al.l ll994'). Emission lines are expected to suffer more screenlike (less grey) 
attenuation than the optical continua, as the emission must first pass through the dense birth 
clouds in which the hot O and B stars that exc ite the nebular emission lines form jCharlot fc 
l2000l : IWild et al.ll2007l : Ida Cunha et al.ll2008l '). An additional complication arises in the case of 
composite galaxies, in which the line emission originating from the narrow line region of the AGN 
may suffer a different level and/or form of extinction compared to nebular emission originating 
in star-forming regions. 

Extinction in the mid-IR is dominated by features caused by the stretching and bending 
modes of amorphous silicate grains. The presenc e of silicate ab sorption can be used to infer the 
geometry of material obscuring the light source ( Le venson et al.ii2007i : iNenkova et al.ll2003 l. 
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Figure 2. Comparing dust obscuration in the IR (10 /im silicate absorption) to dust obscuration 
in the optical (the Balmer decrement). We set B — 2.87 (Case B recombination) for star forming 
and composite galaxies, and B — 3.1 for pure- AGN. Only those galaxies with mean per-pixel 
S/N >4 in the wavelength range 5.4-5.9 /im, and Sl^KS/N >3 in their Ha and H/3 lines are plot- 
ted. The upper axis gives ry, derived from the Balmer decrement assuming the two -component 
Char iot and Fall (2000) dust attenuation curve with ubc = —1.3 and nisM ~ —0.7 l|Wild et al.l 
120071 1. 

In Figure [21 we compare two different measures of dust content in the IR and optical regimes. 
Aside from the two extreme ULIRGs, there is a general trend for the two measures of dust 
content to track one another. Although both axes describe dust content, the trend is slightly 
surprising: any significant silicate absorption applies optically thick clouds (ry > 100), which 
contrasts strongly with the Ty measured here in the optical. Therefore these two quantities 
measure dust in entirely independent regions. It is interesting to note that the two outliers 
from the trend (ULIRGs 15250 and 08572), with extremely strong silicate absorption implying 
that the source of the mid-IR continuum is buried in a smooth, geometrically and optically 
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thick cloud, are composite-AGN in the optical. Therefore some light from the central nucleus 
is escaping to illuminate the narrow-line region, or the AGN is not the dominant light source 
responsible for the mid-IR continuum and the dust is not associated with the AGN torus. 



4. Low-Ionization Emission Lines 

In star-forming galaxies, the dominant excitation mechanism for nebular emission lines is UV 
radiation from young, hot stars. When star formation is the dominant mechanism of excitation. 
Ha is one of the best calibrated measures of instantaneous star formation rate in galaxies. 
However, because the Ha line lies in the optical, to obtain accurate SFR estimates we must 
correct for dust attenuation. In this section we address the question of how accurate the standard 
methods of correction for dust attenuation in the optical are, by comparing the Ha luminosity 
to the luminosity of the lower ionization neon lines in the mid-IR. To first order, the strengths 
of these lines are expected to be strongly correlated, especially in the relatively massive, and 
therefore metal rich, galaxies that make up the majority of our samples. 



10' 



= 10" 



10' 



10" 



m = 2.07 , 



10° 10' 10° 

L(Ha) no dust correction 



10' 



= 10° 



10' 



10° 



m = 0.21 



O K 



■ + • 



10° 10' 10° 10' 10" 

L(Ha) T a X-"' 



10' 



10° 



10' 



10° 



m = 0.73 



X a! 



o ^ .■■ 



10° 



10' 10° 
L(Ha) CFOO + 



10' 



10.0 



1.0 



0.1 



+ SSGSS-HII 1.14 
OSSGSS-comp 1.04 
ASSGSS-AGN 1.89 
A Seyfert 1.78 



i,12112 



A <> 



+ + <*>■ + . 

O-H- 

+ o 



+•+ 



10° 



10' 10° 
L(Ha) CFOO+ 



10' 



Figure 3. The luminosity of Ha in the optical vs. the low ionization Ne lines in the mid-IR. Top 
left: Ha is not corrected for dust attenuation. Top right: Ha is corrected for dust attenuation 
using a single component CFOO model. Bottom left: Ha is corrected for dust attenuation using a 
double component CFOO model. In all three panels, the dashed line shows the best-fit linear fit 
to the SSGSS star-forming galaxies, and the dotted lines show a factor of two above and below. 
Bottom right: The residual from the linear fit in the bottom-left panel. 
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In Figure [3l we plot Hq vs. Ne luminosities for all of the samples, firstly with Ha uncorrected 
for dust attenuation, then corrected with two different dust attenuation laws used in the litera- 
ture. In the top-left panel, where no correction for dust attenuation has been applied to the Ha 
luminosities, we immediately notice that the Balmer-strong galaxies and ULIRGs lie above the 
SSGSS galaxies as expected for their large dust contents. In the top-right panel, the dust law is a 
single component mo del with r oc A~°'^. Th is law provides a good fit to the UV and IR continua 
of starburst galaxies (jCharlot fc "Faill200d . hereafter CFOO), but we see that it does not pro- 
vide a good fit to the nebular emission lines, causing considerable overcorrection of the optical 
emission lines in the dustiest galaxies, and undercorrection in the least dusty Seyferts compared 
to normal SSGSS-Hll galaxies. In the bottom-left panel, we employ the full two-component 
dust model of CFOO, where the second term all ows for the additi o nal, greyer, attenuatio n of the 
stellar birth clouds with an exponent of -1.3 (|Wild et al.ll2007l : IdaCunha et ai]|2008l '). There 
is a noticeable improvement in the scatter. In each panel, the dashed line shows a linear least- 
squares fit to the SSGSS-Hll galax;ies, with the dotted lines indicating a factor of two above and 
below the best-fit. In the bottom-right we plot the residual away from this best-fit line using the 
two-component CFOO model. The scatter is around a factor of two, even for the extreme dusty 
Balmer-strong galaxies. Even the ULIRGs do not show a dramatic underprediction of their Ha 
luminosities, despite their optical Balmer decrements indicating as much as six magnitudes of 
extinction and their silicate absorption indicating hundreds of magnitudes. We note a slight 
mean offset for the AGN, in the sense that they have stronger Nell-I-Nelll emission for their Ha 
luminosities. This could be due to the harder ionizing radiation field of the AGN, or that the 
true dust attenuation law is slightly different in the narrow line region of the AGN. Although 
the effect is not large, further study may be able to untangle these two effects. 



5. High-Ionization Emission Lines 

Given the success of the 2 component dust law for correcting Ha for dust attenuation in 
even the dustiest galaxies in the local Universe, we now turn to the question of using [O ill] 
luminosity as a proxy for black hole accretion rate. Although originating from the narrow line 
region, far removed spatially from the central nucleus with its dense obscuring torus, the [O ill] 
emission must still suffer extinction from the host galaxy, and possibly also from dust within 
the narrow line region. This adds an additional degree of complexity to the problem: the Ha 
and H/3 emission lines with which we calculate the magnitude of the dust attenuation in the 
optical, may originate at least in part from star-forming regions completely spatially separated 
from the narrow-line region in which we are interested. However, given the small residual offset 
between Ne and Ha emission presented in the previous section, even for the strongest AGN in 
the SDSS (the Seyfert sample), this effect cannot be large. 

In Figure |3J we plot individual detections of JO iv] A25.9 versus [O ill] A5007, using the 
Balmer decrement and two component CFOO model to correct [O ill] for dust attenuation. The 
mid-IR [O iv] line lies too close to an Fe ll line for deconvolution in the low resolution Spttzer 
spectra. Therefore we only plot the the Seyfert, ULIRG and dusty Balmer-strong galaxy samples, 
for which we have high resolution spectra. Of the galaxies in this figure, two of the Balmer-strong 
galaxies and both of the ULIRGs are optically classified as pure-AGN, and the majority of the 
[O III] is therefore likely to originate from the AGN rather than from star formation. For the 
third Balmer-strong galaxy, a small correction has been applied to account for the star formation 
contamination of the [O ill] line. At higher ionization, the [O iv] line is expected to suffer less 
contamination from star formation. 

Once again, we find that where there are individual detections of both [O iv] and [O ill] , 
the line luminosities are strongly correlated after appropriate correction for dust attenuation of 
[O III] , even in extremely dusty galaxies. Although further work is required to investigate those 
galaxies where [O iv] is not detected in the Spitzer spectra, this preliminary result suggests that 
dust attenuation is not a significant cause for concern in the use of optical spectra for deriving 
black hole accretion rates. 
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Figure 4. Dust-corrected L([Olll]) vs. L([0iv]) for the Seyfert sample (filled triangles), the 
three Balmer-strong galaxies with [O iv] detections (red circles) and the two ULIRGs with both 
[O iv] detections and SDSS spectra. 



6. Discussion 

In this preliminary investigation comparing the mid-IR and optical emission lines of three 
samples of unusual galaxies in the local Universe, and a sample of more ordinary galaxies, we 
find little evidence that the optical wavelength regime, used alone, would provide a misleading 
picture in terms of either star formation rates or black hole accretion rates. Each of the three 
samples is extreme in some way: the Seyferts are the strongest [O ill] emitters in the local 
Universe, the Balmer-strong galaxies are some of the dustiest AGN in the SDSS, the ULIRGs 
are unusually bright in the IR, implying substantial dust contents. We have shown how each of 
these samples fits into the global galaxy population, both in the mid-IR and in the optical. 

It is perhaps surprising that, given the geometrically complicated dust distributions of some 
galaxies, and the spatial separation between the narrow line region and the star-forming re- 
gions in composite- AGN, that nebular emission lines in the optical wavelength regime can be 
adequately corrected for the effects of dust attenuation. With further investigation of these sam- 
ples, we hope to understand more about the limitations, or otherwise, of the optical wavelength 
regime for studying unusual galaxies. 
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